Functional annotation signatures of disease susceptibility loci improve SNP association analysis.

BACKGROUND: Genetic association studies are conducted to discover genetic loci that contribute to an inherited trait, identify the variants behind these associations and ascertain their functional role in determining the phenotype. To date, functional annotations of the genetic variants have rarely played more than an indirect role in assessing evidence for association. Here, we demonstrate how these data can be systematically integrated into an association study's analysis plan. RESULTS: We developed a Bayesian statistical model for the prior probability of phenotype-genotype association that incorporates data from past association studies and publicly available functional annotation data regarding the susceptibility variants under study. The model takes the form of a binary regression of association status on a set of annotation variables whose coefficients were estimated through an analysis of associated SNPs in the GWAS Catalog (GC). The functional predictors examined included measures that have been demonstrated to correlate with the association status of SNPs in the GC and some whose utility in this regard is speculative: summaries of the UCSC Human Genome Browser ENCODE super-track data, dbSNP function class, sequence conservation summaries, proximity to genomic variants in the Database of Genomic Variants and known regulatory elements in the Open Regulatory Annotation database, PolyPhen-2 probabilities and RegulomeDB categories. Because we expected that only a fraction of the annotations would contribute to predicting association, we employed a penalized likelihood method to reduce the impact of non-informative predictors and evaluated the model's ability to predict GC SNPs not used to construct the model. We show that the functional data alone are predictive of a SNP's presence in the GC. Further, using data from a genome-wide study of ovarian cancer, we demonstrate that their use as prior data when testing for association is practical at the genome-wide scale and improves power to detect associations. CONCLUSIONS: We show how diverse functional annotations can be efficiently combined to create 'functional signatures' that predict the a priori odds of a variant's association to a trait and how these signatures can be integrated into a standard genome-wide-scale association analysis, resulting in improved power to detect truly associated variants.

Role of mast cells in inflammatory bowel disease and inflammation-associated colorectal neoplasia in IL-10-deficient mice.

BACKGROUND: Inflammatory bowel disease (IBD) is hypothesized to result from stimulation of immune responses against resident intestinal bacteria within a genetically susceptible host. Mast cells may play a critical role in IBD pathogenesis, since they are typically located just beneath the intestinal mucosal barrier and can be activated by bacterial antigens. METHODOLOGY/PRINCIPAL FINDINGS: This study investigated effects of mast cells on inflammation and associated neoplasia in IBD-susceptible interleukin (IL)-10-deficient mice with and without mast cells. IL-10-deficient mast cells produced more pro-inflammatory cytokines in vitro both constitutively and when triggered, compared with wild type mast cells. However despite this enhanced in vitro response, mast cell-sufficient Il10(-/-) mice actually had decreased cecal expression of tumor necrosis factor (TNF) and interferon (IFN)-gamma mRNA, suggesting that mast cells regulate inflammation in vivo. Mast cell deficiency predisposed Il10(-/-) mice to the development of spontaneous colitis and resulted in increased intestinal permeability in vivo that preceded the development of colon inflammation. However, mast cell deficiency did not affect the severity of IBD triggered by non-steroidal anti-inflammatory agents (NSAID) exposure or helicobacter infection that also affect intestinal permeability. CONCLUSIONS/SIGNIFICANCE: Mast cells thus appear to have a primarily protective role within the colonic microenvironment by enhancing the efficacy of the mucosal barrier. In addition, although mast cells were previously implicated in progression of sporadic colon cancers, mast cells did not affect the incidence or severity of colonic neoplasia in this inflammation-associated model.

Dusp3 and Psme3 are associated with murine susceptibility to Staphylococcus aureus infection and human sepsis.

Using A/J mice, which are susceptible to Staphylococcus aureus, we sought to identify genetic determinants of susceptibility to S. aureus, and evaluate their function with regard to S. aureus infection. One QTL region on chromosome 11 containing 422 genes was found to be significantly associated with susceptibility to S. aureus infection. Of these 422 genes, whole genome transcription profiling identified five genes (Dcaf7, Dusp3, Fam134c, Psme3, and Slc4a1) that were significantly differentially expressed in a) S. aureus -infected susceptible (A/J) vs. resistant (C57BL/6J) mice and b) humans with S. aureus blood stream infection vs. healthy subjects. Three of these genes (Dcaf7, Dusp3, and Psme3) were down-regulated in susceptible vs. resistant mice at both pre- and post-infection time points by qPCR. siRNA-mediated knockdown of Dusp3 and Psme3 induced significant increases of cytokine production in S. aureus-challenged RAW264.7 macrophages and bone marrow derived macrophages (BMDMs) through enhancing NF-κB signaling activity. Similar increases in cytokine production and NF-κB activity were also seen in BMDMs from CSS11 (C57BL/6J background with chromosome 11 from A/J), but not C57BL/6J. These findings suggest that Dusp3 and Psme3 contribute to S. aureus infection susceptibility in A/J mice and play a role in human S. aureus infection.

Lysosomal Disorders Drive Susceptibility to Tuberculosis by Compromising Macrophage Migration.

A zebrafish genetic screen for determinants of susceptibility to Mycobacterium marinum identified a hypersusceptible mutant deficient in lysosomal cysteine cathepsins that manifests hallmarks of human lysosomal storage diseases. Under homeostatic conditions, mutant macrophages accumulate undigested lysosomal material, which disrupts endocytic recycling and impairs their migration to, and thus engulfment of, dying cells. This causes a buildup of unengulfed cell debris. During mycobacterial infection, macrophages with lysosomal storage cannot migrate toward infected macrophages undergoing apoptosis in the tuberculous granuloma. The unengulfed apoptotic macrophages undergo secondary necrosis, causing granuloma breakdown and increased mycobacterial growth. Macrophage lysosomal storage similarly impairs migration to newly infecting mycobacteria. This phenotype is recapitulated in human smokers, who are at increased risk for tuberculosis. A majority of their alveolar macrophages exhibit lysosomal accumulations of tobacco smoke particulates and do not migrate to Mycobacterium tuberculosis. The incapacitation of highly microbicidal first-responding macrophages may contribute to smokers' susceptibility to tuberculosis.

Animal models of soft-tissue sarcoma.

Soft-tissue sarcomas (STSs) are rare mesenchymal tumors that arise from muscle, fat and connective tissue. Currently, over 75 subtypes of STS are recognized. The rarity and heterogeneity of patient samples complicate clinical investigations into sarcoma biology. Model organisms might provide traction to our understanding and treatment of the disease. Over the past 10 years, many successful animal models of STS have been developed, primarily genetically engineered mice and zebrafish. These models are useful for studying the relevant oncogenes, signaling pathways and other cell changes involved in generating STSs. Recently, these model systems have become preclinical platforms in which to evaluate new drugs and treatment regimens. Thus, animal models are useful surrogates for understanding STS disease susceptibility and pathogenesis as well as for testing potential therapeutic strategies.

A longitudinal study of epigenetic variation in twins.

DNA methylation is a key epigenetic mechanism involved in the developmental regulation of gene expression. Alterations in DNA methylation are established contributors to inter-individual phenotypic variation and have been associated with disease susceptibility. The degree to which changes in loci-specific DNA methylation are under the influence of heritable and environmental factors is largely unknown. In this study, we quantitatively measured DNA methylation across the promoter regions of the dopamine receptor 4 gene (DRD4), the serotonin transporter gene (SLC6A4/SERT) and the X-linked monoamine oxidase A gene (MAOA) using DNA sampled at both ages 5 and 10 years in 46 MZ twin-pairs and 45 DZ twin-pairs (total n=182). Our data suggest that DNA methylation differences are apparent already in early childhood, even between genetically identical individuals, and that individual differences in methylation are not stable over time. Our longitudinal-developmental study suggests that environmental influences are important factors accounting for interindividual DNA methylation differences, and that these influences differ across the genome. The observation of dynamic changes in DNA methylation over time highlights the importance of longitudinal research designs for epigenetic research.

High-throughput isolation and mapping of C. elegans mutants susceptible to pathogen infection.

We present a novel strategy that uses high-throughput methods of isolating and mapping C. elegans mutants susceptible to pathogen infection. We show that C. elegans mutants that exhibit an enhanced pathogen accumulation (epa) phenotype can be rapidly identified and isolated using a sorting system that allows automation of the analysis, sorting, and dispensing of C. elegans by measuring fluorescent bacteria inside the animals. Furthermore, we validate the use of Amplifluor as a new single nucleotide polymorphism (SNP) mapping technique in C. elegans. We show that a set of 9 SNPs allows the linkage of C. elegans mutants to a 5-8 megabase sub-chromosomal region.

WHO guidelines for antimicrobial treatment in children admitted to hospital in an area of intense Plasmodium falciparum transmission: prospective study.

OBJECTIVES: To assess the performance of WHO's "Guidelines for care at the first-referral level in developing countries" in an area of intense malaria transmission and identify bacterial infections in children with and without malaria. DESIGN: Prospective study. SETTING: District hospital in Muheza, northeast Tanzania. PARTICIPANTS: Children aged 2 months to 13 years admitted to hospital for febrile illness. MAIN OUTCOME MEASURES: Sensitivity and specificity of WHO guidelines in diagnosing invasive bacterial disease; susceptibility of isolated organisms to recommended antimicrobials. RESULTS: Over one year, 3639 children were enrolled and 184 (5.1%) died; 2195 (60.3%) were blood slide positive for Plasmodium falciparum, 341 (9.4%) had invasive bacterial disease, and 142 (3.9%) were seropositive for HIV. The prevalence of invasive bacterial disease was lower in slide positive children (100/2195, 4.6%) than in slide negative children (241/1444, 16.7%). Non-typhi Salmonella was the most frequently isolated organism (52/100 (52%) of organisms in slide positive children and 108/241 (45%) in slide negative children). Mortality among children with invasive bacterial disease was significantly higher (58/341, 17%) than in children without invasive bacterial disease (126/3298, 3.8%) (P<0.001), and this was true regardless of the presence of P falciparum parasitaemia. The sensitivity and specificity of WHO criteria in identifying invasive bacterial disease in slide positive children were 60.0% (95% confidence interval 58.0% to 62.1%) and 53.5% (51.4% to 55.6%), compared with 70.5% (68.2% to 72.9%) and 48.1% (45.6% to 50.7%) in slide negative children. In children with WHO criteria for invasive bacterial disease, only 99/211(47%) of isolated organisms were susceptible to the first recommended antimicrobial agent. CONCLUSIONS: In an area exposed to high transmission of malaria, current WHO guidelines failed to identify almost a third of children with invasive bacterial disease, and more than half of the organisms isolated were not susceptible to currently recommended antimicrobials. Improved diagnosis and treatment of invasive bacterial disease are needed to reduce childhood mortality.

Dissecting the Functional Impacts of Non-Coding Genetic Variation

A large proportion of the variation in traits between individuals can be attributed to variation in the nucleotide sequence of the genome. The most commonly studied traits in human genetics are related to disease and disease susceptibility. Although scientists have identified genetic causes for over 4,000 monogenic diseases, the underlying mechanisms of many highly prevalent multifactorial inheritance disorders such as diabetes, obesity, and cardiovascular disease remain largely unknown. Identifying genetic mechanisms for complex traits has been challenging because most of the variants are located outside of protein-coding regions, and determining the effects of such non-coding variants remains difficult. In this dissertation, I evaluate the hypothesis that such non-coding variants contribute to human traits and diseases by altering the regulation of genes rather than the sequence of those genes. I will specifically focus on studies to determine the functional impacts of genetic variation associated with two related complex traits: gestational hyperglycemia and fetal adiposity. At the genomic locus associated with maternal hyperglycemia, we found that genetic variation in regulatory elements altered the expression of the HKDC1 gene. Furthermore, we demonstrated that HKDC1 phosphorylates glucose in vitro and in vivo, thus demonstrating that HKDC1 is a fifth human hexokinase gene. At the fetal-adiposity associated locus, we identified variants that likely alter VEPH1 expression in preadipocytes during differentiation. To make such studies of regulatory variation high-throughput and routine, we developed POP-STARR, a novel high throughput reporter assay that can empirically measure the effects of regulatory variants directly from patient DNA. By combining targeted genome capture technologies with STARR-seq, we assayed thousands of haplotypes from 760 individuals in a single experiment. We subsequently used POP-STARR to identify three key features of regulatory variants: that regulatory variants typically have weak effects on gene expression; that the effects of regulatory variants are often coordinated with respect to disease-risk, suggesting a general mechanism by which the weak effects can together have phenotypic impact; and that nucleotide transversions have larger impacts on enhancer activity than transitions. Together, the findings presented here demonstrate successful strategies for determining the regulatory mechanisms underlying genetic associations with human traits and diseases, and value of doing so for driving novel biological discovery.

A genome-wide association study identifies susceptibility variants for type 2 diabetes in Han Chinese.

To investigate the underlying mechanisms of T2D pathogenesis, we looked for diabetes susceptibility genes that increase the risk of type 2 diabetes (T2D) in a Han Chinese population. A two-stage genome-wide association (GWA) study was conducted, in which 995 patients and 894 controls were genotyped using the Illumina HumanHap550-Duo BeadChip for the first genome scan stage. This was further replicated in 1,803 patients and 1,473 controls in stage 2. We found two loci not previously associated with diabetes susceptibility in and around the genes protein tyrosine phosphatase receptor type D (PTPRD) (P = 8.54x10(-10); odds ratio [OR] = 1.57; 95% confidence interval [CI] = 1.36-1.82), and serine racemase (SRR) (P = 3.06x10(-9); OR = 1.28; 95% CI = 1.18-1.39). We also confirmed that variants in KCNQ1 were associated with T2D risk, with the strongest signal at rs2237895 (P = 9.65x10(-10); OR = 1.29, 95% CI = 1.19-1.40). By identifying two novel genetic susceptibility loci in a Han Chinese population and confirming the involvement of KCNQ1, which was previously reported to be associated with T2D in Japanese and European descent populations, our results may lead to a better understanding of differences in the molecular pathogenesis of T2D among various populations.

Two genes on A/J chromosome 18 are associated with susceptibility to Staphylococcus aureus infection by combined microarray and QTL analyses.

Although it has recently been shown that A/J mice are highly susceptible to Staphylococcus aureus sepsis as compared to C57BL/6J, the specific genes responsible for this differential phenotype are unknown. Using chromosome substitution strains (CSS), we found that loci on chromosomes 8, 11, and 18 influence susceptibility to S. aureus sepsis in A/J mice. We then used two candidate gene selection strategies to identify genes on these three chromosomes associated with S. aureus susceptibility, and targeted genes identified by both gene selection strategies. First, we used whole genome transcription profiling to identify 191 (56 on chr. 8, 100 on chr. 11, and 35 on chr. 18) genes on our three chromosomes of interest that are differentially expressed between S. aureus-infected A/J and C57BL/6J. Second, we identified two significant quantitative trait loci (QTL) for survival post-infection on chr. 18 using N(2) backcross mice (F(1) [C18A]xC57BL/6J). Ten genes on chr. 18 (March3, Cep120, Chmp1b, Dcp2, Dtwd2, Isoc1, Lman1, Spire1, Tnfaip8, and Seh1l) mapped to the two significant QTL regions and were also identified by the expression array selection strategy. Using real-time PCR, 6 of these 10 genes (Chmp1b, Dtwd2, Isoc1, Lman1, Tnfaip8, and Seh1l) showed significantly different expression levels between S. aureus-infected A/J and C57BL/6J. For two (Tnfaip8 and Seh1l) of these 6 genes, siRNA-mediated knockdown of gene expression in S. aureus-challenged RAW264.7 macrophages induced significant changes in the cytokine response (IL-1 beta and GM-CSF) compared to negative controls. These cytokine response changes were consistent with those seen in S. aureus-challenged peritoneal macrophages from CSS 18 mice (which contain A/J chromosome 18 but are otherwise C57BL/6J), but not C57BL/6J mice. These findings suggest that two genes, Tnfaip8 and Seh1l, may contribute to susceptibility to S. aureus in A/J mice, and represent promising candidates for human genetic susceptibility studies.

Genome-wide scan of copy number variation in late-onset Alzheimer's disease.

Alzheimer's disease is a complex and progressive neurodegenerative disease leading to loss of memory, cognitive impairment, and ultimately death. To date, six large-scale genome-wide association studies have been conducted to identify SNPs that influence disease predisposition. These studies have confirmed the well-known APOE epsilon4 risk allele, identified a novel variant that influences disease risk within the APOE epsilon4 population, found a SNP that modifies the age of disease onset, as well as reported the first sex-linked susceptibility variant. Here we report a genome-wide scan of Alzheimer's disease in a set of 331 cases and 368 controls, extending analyses for the first time to include assessments of copy number variation. In this analysis, no new SNPs show genome-wide significance. We also screened for effects of copy number variation, and while nothing was significant, a duplication in CHRNA7 appears interesting enough to warrant further investigation.

Overexpression of the cardiac beta(2)-adrenergic receptor and expression of a beta-adrenergic receptor kinase-1 (betaARK1) inhibitor both increase myocardial contractility but have differential effects on susceptibility to ischemic injury.

Cardiac beta(2)-adrenergic receptor (beta(2)AR) overexpression is a potential contractile therapy for heart failure. Cardiac contractility was elevated in mice overexpressing beta(2)ARs (TG4s) with no adverse effects under normal conditions. To assess the consequences of beta(2)AR overexpression during ischemia, perfused hearts from TG4 and wild-type mice were subjected to 20-minute ischemia and 40-minute reperfusion. During ischemia, ATP and pH fell lower in TG4 hearts than wild type. Ischemic injury was greater in TG4 hearts, as indicated by lower postischemic recoveries of contractile function, ATP, and phosphocreatine. Because beta(2)ARs, unlike beta(1)ARs, couple to G(i) as well as G(s), we pretreated mice with the G(i) inhibitor pertussis toxin (PTX). PTX treatment increased basal contractility in TG4 hearts and abolished the contractile resistance to isoproterenol. During ischemia, ATP fell lower in TG4+PTX than in TG4 hearts. Recoveries of contractile function and ATP were lower in TG4+PTX than in TG4 hearts. We also studied mice that overexpressed either betaARK1 (TGbetaARK1) or a betaARK1 inhibitor (TGbetaARKct). Recoveries of function, ATP, and phosphocreatine were higher in TGbetaARK1 hearts than in wild-type hearts. Despite basal contractility being elevated in TGbetaARKct hearts to the same level as that of TG4s, ischemic injury was not increased. In summary, beta(2)AR overexpression increased ischemic injury, whereas betaARK1 overexpression was protective. Ischemic injury in the beta(2)AR overexpressors was exacerbated by PTX treatment, implying that it was G(s) not G(i) activity that enhanced injury. Unlike beta(2)AR overexpression, basal contractility was increased by betaARK1 inhibitor expression without increasing ischemic injury, thus implicating a safer potential therapy for heart failure.

Polymorphisms at the microRNA binding-site of the stem cell marker gene CD133 modify susceptibility to and survival of gastric cancer.

CD133 is one of the most common stem cell markers, and functional single nucleotide polymorphisms (SNPs) of CD133 may modulate its gene functions and thus cancer risk and patient survival. We hypothesized that potentially functional CD133 SNPs are associated with gastric cancer (GC) risk and survival. To test this hypothesis, we conducted a case-control study of 371 GC patients and 313 cancer-free controls frequency-matched by age, sex, and ethnicity. We genotyped four selected, potentially functional CD133 SNPs (rs2240688A>C, rs7686732C>G, rs10022537T>A, and rs3130C>T) and used logistic regression analysis for associations of these SNPs with GC risk and Cox hazards regression analysis for survival. We found that compared with the miRNA binding site rs2240688 AA genotype, AC + CC genotypes were associated with significantly increased GC risk (adjusted OR = 1.52, 95% CI = 1.09-2.13); for another miRNA binding site rs3130C>T SNP, the TT genotype was associated with significantly reduced GC risk (adjusted OR = 0.68, 95% CI = 0.48-0.97), compared with CC + CT genotypes. In all patients, the risk rs3130 TT variant genotype was significantly associated with overall survival (OS) (adjusted P(trend) = 0.016 and 0.007 under additive and recessive models, respectively). These findings suggest that these two CD133 miRNA binding site variants, rs2240688 and rs3130, may be potential biomarkers for genetic susceptibility to GC and possible predictors for survival in GC patients but require further validation by larger studies.

Probing Tissue Microstructure Using Susceptibility Contrast Magnetic Resonance Imaging

Magnetic resonance imaging is a research and clinical tool that has been applied in a wide variety of sciences. One area of magnetic resonance imaging that has exhibited terrific promise and growth in the past decade is magnetic susceptibility imaging. Imaging tissue susceptibility provides insight into the microstructural organization and chemical properties of biological tissues, but this image contrast is not well understood. The purpose of this work is to develop effective approaches to image, assess, and model the mechanisms that generate both isotropic and anisotropic magnetic susceptibility contrast in biological tissues, including myocardium and central nervous system white matter.

This document contains the first report of MRI-measured susceptibility anisotropy in myocardium. Intact mouse heart specimens were scanned using MRI at 9.4 T to ascertain both the magnetic susceptibility and myofiber orientation of the tissue. The susceptibility anisotropy of myocardium was observed and measured by relating the apparent tissue susceptibility as a function of the myofiber angle with respect to the applied magnetic field. A multi-filament model of myocardial tissue revealed that the diamagnetically anisotropy α-helix peptide bonds in myofilament proteins are capable of producing bulk susceptibility anisotropy on a scale measurable by MRI, and are potentially the chief sources of the experimentally observed anisotropy.

The growing use of paramagnetic contrast agents in magnetic susceptibility imaging motivated a series of investigations regarding the effect of these exogenous agents on susceptibility imaging in the brain, heart, and kidney. In each of these organs, gadolinium increases susceptibility contrast and anisotropy, though the enhancements depend on the tissue type, compartmentalization of contrast agent, and complex multi-pool relaxation. In the brain, the introduction of paramagnetic contrast agents actually makes white matter tissue regions appear more diamagnetic relative to the reference susceptibility. Gadolinium-enhanced MRI yields tensor-valued susceptibility images with eigenvectors that more accurately reflect the underlying tissue orientation.

Despite the boost gadolinium provides, tensor-valued susceptibility image reconstruction is prone to image artifacts. A novel algorithm was developed to mitigate these artifacts by incorporating orientation-dependent tissue relaxation information into susceptibility tensor estimation. The technique was verified using a numerical phantom simulation, and improves susceptibility-based tractography in the brain, kidney, and heart. This work represents the first successful application of susceptibility-based tractography to a whole, intact heart.

The knowledge and tools developed throughout the course of this research were then applied to studying mouse models of Alzheimer’s disease in vivo, and studying hypertrophic human myocardium specimens ex vivo. Though a preliminary study using contrast-enhanced quantitative susceptibility mapping has revealed diamagnetic amyloid plaques associated with Alzheimer’s disease in the mouse brain ex vivo, non-contrast susceptibility imaging was unable to precisely identify these plaques in vivo. Susceptibility tensor imaging of human myocardium specimens at 9.4 T shows that susceptibility anisotropy is larger and mean susceptibility is more diamagnetic in hypertrophic tissue than in normal tissue. These findings support the hypothesis that myofilament proteins are a source of susceptibility contrast and anisotropy in myocardium. This collection of preclinical studies provides new tools and context for analyzing tissue structure, chemistry, and health in a variety of organs throughout the body.